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O '. Abstract 
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Light dark matter (DM) models with isospin violation (ISV) and large DM-nucleon spin- 
independent cross section (T„ may provide a way to understand the confusing DM direct detection 
experimental results. Combining with stringent astrophysical and collider constraints, we can fur- 
ther deduce the DM properties. In light of general operator analyses, we systematically study the 
ISV operators with fn/ fp < 0, and show that the required ISV fn/ fp — —0.7 must arise from the 
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o ■ DM and first-family quark couplings. Then we discuss three kinds of the ISV models classified 

(— I ' by the mediators: a light Z' gauge boson in the extra U{\)x model, a (approximate) spectator 

Higgs doublet, and color triplets. In particular, the spectator Higgs doublet model can explain the 
^ ' Tevatron CDF W + jj anomaly simultaneously. In addition, although the U{l)x gauge boson Z' 

^vq I which has kinetic mixing with U{1)y gauge boson only generates /„ = 0, we can combine it with 

the conventional Higgs Yukawa couplings to achieve the proper ISV. Especially, most of our models 
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t~^ ■ can address the recent GoGeNT annual modulation as well as all the other DM direct detection 

o ■ 

experiments. As a concrete example, we propose the U{l)x model where the U{l)x charged light 
sneutrino is the inelastic DM (iDM), and dominantly decays to the light dark states such as Z' 

ij {Mz' <1 GeV). With ISV it is consistent with all the DM direct detection experiments and satisfies 

rS 

H ' all the other constraints. 
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I. INTRODUCTION AND MOTIVATION 

The way that dark matter (DM) interacts with nucleons is a puzzle by virtue of the 
absence of confirmative experiments results. The recent possible progress made on direct 
detections may shed light on it. On the one hand, the events from the DAM A/LIBRA [1] 
and CoGeNT [2] experiments imply a light DM (LDM) with mass around 8 GeV and a rather 
large spin-independent (SI) cross section with nucleon o"„ ~ 2 x 10""^ pb and (T„ ~ 5 x 10~^ 
pb, respectively. On the other hand, the null results from XENON [4] and CDMS [3] 
experiments challenge the CoGeNT/DAMA results. Moreover, for the CoGeNT/DAMA 
favored DM there may be some expected signals (but not found) from astrophysics or collider 
experiments. So how to reconcile these results in a natural way may guide us to identify 
some of DM properties. 

DM-nucleus interactions with isospin-violation (ISV) may provide an interesting way to 
reconcile various direct detection experimental results [5]. ISV is not a novel phenomena [6, 
7], which arises when the DM interactions with proton and neutron have different strengths, 
fp ^ In- If fn/ fp < the amplitudes between the DM-proton and DM- neutron destructively 
interfere, leading to a cancellation in the DM-nucleus amplitude. The degree of cancellation 
varies with the concrete nucleus used in the detector target. In this way, the strongest 
constraints from XENON are substantially weakened by taking fn/ fp = —0.7, especially, 
the CoGeNT/DAMA regions may overlap [5]. But this region is excluded by the CDMS-Ge 
experiment which use the same atom as the CoGeNT experiment. 

However, if one only considers annual modulation (the observed CoGeNT annual mod- 
ulation with significance of 2.8a [8]), inelastic DM (iDM) scenario [9] is able to enhance 
modulation and reduce tension between CDMS-Ge and CoGeNT. The Ref. [10] found that 
taking fn/ fp = —0.7 and the upper-limit quenching factor Qj^a = 0.43, one can address all 
the confusing experimental results consistently via an iDM ~10 GeV with mass splitting 
5^15 keV. On the other hand, right now the direct detection experiments may have reached 
the level to distinguish /„ and fp [1 1] . Once determined, they definitely convey information 
on the DM-quark interactions. In a word, it is worthy of studying the LDM models with 
ISV systematically. 

Inspired by the aforementioned facts, in this work we make rather thorough analyses on 
the LDM models with ISV, focussing on the scalar and fermionic DMs. We outline the 
generic characters for the LDM models: 

• Proper ISV fn/fp ~ —0.7. How to get this ISV is highly non-trivial. In fact, we shall 
show that the conventional mediators, such as the Higgs/Z/squarks fail to accommo- 
date this value at least in the case of a single dominated mediator. So, we need to 
investigate the new mediators beyond them. 



• Large DM-nucleon SI scattering cross section o"„ ~ 10~^ pb. In the ISV models, 
the nucleus amplitude is reduced by destructively interference. So cr„ is required to 
be about 2 orders larger than the conventional scenario. However, such large cross 
section will bring tension (or even be excluded) with some astrophysical or collider 
constraints. 

• Right DM relic density QDuh"^ — 0.11. We may relax this constraint by allowing 
possible new annihilation channels beyond the Standard Model (SM) fermions as final 
states. In actual model building, LDM usually can annihilate to light dark sector 
states, which will be very helpful sometime. 

• The possible minimal models which may explain the other new physics as well. Espe- 
cially, we try to account for the Tavetron physics such as the CDF W+2jets anoma- 
lies [12]. 

In this paper, we systematically study the ISV operators with fn/ fp < via the effective 
operaror analyses. In particular, we show that the required ISV fn/ fp — —0.7 must arise 
from the couplings between the DM and first-family quarks. Moreover, we propose three 
kinds of the ISV models classified by the mediators: a light Z' gauge boson in the extra 
U{l)x model, a (approximate) spectator Higgs doublet, and color triplets. Interestingly, the 
spectator Higgs doublet model can explain the Tevatron W + jj anomaly simultaneously. In 
addition, combining the U{l)x gauge boson Z' which has kinetic mixing with U{1)y gauge 
boson with the conventional Higgs Yukawa couplings, we also obtain the proper ISV. We 
emphasize that most of our models can address the recent GoGeNT annual modulation as 
well as all the other DM direct detection experiments. To be concrete, we propose the U{l)x 
model where the U{l)x charged light sneutrino is the inelastic DM (iDM), and dominantly 
decays to the light dark states such as Z' {Mz' <1 GeV). With ISV it is not only consistent 
with all the DM direct detection experiments but also satisfies all the other constraints. 

The paper is organized as follows. In Section II, we make some general effective operator 
analyses for the dark matter models with ISV. In Section HI, we present three kinds of the 
ISV models by considering the different mediators. In Section IV, a concrete model with 
sneutrino iDM is discussed. Conclusion and discussions are given in Section V, and finally 
some useful formulas are collected in the Appendices. 

II. GENERIC OPERATOR ANALYSES 

Let us address the setup for generic operator analyses which simplify the discussions. 
The ISV mainly comes from the different U{\)x charges for ur and d^ quarks, or from the 



different Yukawa couplings between the DM and the first-family quarks. 

A. Effective ISV Operators for the CoGeNT/DAMA Experiments 

General anatomy of the interactions between the DM and SM fields involves a large class 
of effective operators Or/uf^f with F = 1,7^,75... where / is the SM fermion, and Cdm 
denote the corresponding DM bilinear operators. However, they are greatly reduced if we are 
just interested in the operators that are relevant to the CoGeNT/DAMA experiments. And 
we can further recover their SU{3)c x U{l)EM—^y completion by specifying the possible 
mediators that connect the DM and the SM fermions, please see the Appendix A for details. 

At first, we pick out the operators generating SI cross sections, and only three of them 
are not non-relativistic (NR) suppressed [13, 14]: 

qq, qii^q, {G%y. (1) 

In a concrete model, only a subset of them maybe generated simultaneously. Obviously, the 
gluon operators can not generate ISV and then are dropped in our discussions. Converting 
to the operators on the DM-SM particle interactions, we obtain the scalar and vector inter- 
actions, which are respectively given in Eqs. (2) and (3). In this paper, we use a^ and bq to 
describe the operator coefficients for the scalar and vector type interactions, respectively. 

aq.XXMh aq^4>^ (pqiqi, (2) 

bqiXl'^XQili^Qi, Krd''(l)qai,qi, (3) 

where x is a fermionic DM while is a scalar DM. If x is Majorana or is real, then the cor- 
responding vector interactions vanishes automatically. Operators xi^XQQ a^^id Xl^l^XQl^Q. 
also generate SI cross section, nevertheless NR suppressed. In our notation, q = u,d denotes 
the up- and down-type quarks respectively, with i the family index (all fields are written in 
the mass eigenstates). We will show that only the first-family quarks are crucial to produce 
ISV. To suppress the potential large flavor violation, the quark bilinears are supposed to 
be diagonal in the flavor space. But in some cases, we will find that this assumption is not 
hold so naively. The operator coefficients a^ and bg are assumed to be suppressed by high 
threshold scale (e.g., mass scale of mediator) A ^ Mdm unless we specify. 

Next, there are other operators that do not contribute to the leading order SI cross section 
in the NR limit, nevertheless usually generated together with the other operators in Eq. (2) 
and Eq. (3) in a UV completed theory. Interestingly, they are important to determine the 
relic density of DM. The indirect detection signals and the annihilation rates from these 
operators are collected in the Appendix A. To be concrete, we present these operators in 



the following 

xi^xfi^f, xxh^f, xi^xff, 

xi^if^xh^Yf, xl^.xh^Yf, Xl^l^^xfYf, 

I0P/75/, 0^9^0/7^/- • (4) 

For the LDM we only have to consider the SM states lighter than Mdm which the LDM 
can annihilate into. So Eqs. (2)- (4) give a general effective operator description on the DM 
models inspired by the CoGeNT/DAMA experiments. Having established the setup, in the 
following we shall investigate the ISV from scalar and vector interactions. 

First let us consider scalar interactions given in Eq. (2). To construct the DM-nucleon 
effective operators from microscopic DM-quark interactions, we must calculate the quark 
bilinear matrix elements in the nucleon states. In the case of scalar interaction, they are 
given by 

mq{n\qq\n) = mnf^\ (5) 

where n denotes either the proton or neutron, and nig the quark mass. For the light quarks, 
the form factors f^^ = 0.020 ± 0.004, f^^ = 0.014 ± 0.003, f^^J = 0.026 ± 0.005, 4^^ = 
0.036 ± 0.008, and /^^'"^ = 0.118 ± 0.062 [15]. The heavy quarks q = c,b,t contribute to the 
nucleon mass through the triangle diagram [16] as follows 




mg{n\qq\n) = —m^ 1 - V fr\= T^rrinff;, (6) 



where Jtq — 0.84 and 0.83 for the proton and neutron respectively. Through the above two 
equations, we get the DM-nucleon effective couplings 



/„ OC a„ = TTlr. 
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c,b,t ^. 



g=quarks 



where the dimensionless quantities B'^ = /^"^m„/mg, encoding the ISV in the nucleus itself. 
We use the quark masses: m„ = 0.002 GeV, m^ = 0.005 GeV, m^ = 0.095 GeV, rric = 1.25 
GeV, rrih = 4.2 GeV, rrit = 172.3 GeV. Then we get 

5^^9.3, B:^Q.5, 5^^5.1, 5,-^7.1, (8) 

5P'"^1.2, BP'"^0.05, fi^'"^ 0.015, fif'"- 0.00035. (9) 

By the way, in Eq. (7) /„ = «„ for a fermionic DM. Whereas for s scalar DM a^ has dimension 
+1, and if the DM-nucleus SI scattering cross section written in the form of Eq. (Bl), then 
actually fn = a„/2MDM- Hereafter we will absorb the 1/2Mdm into a^ for the scalar DM. 



Remarkably, from Eqs. (7)- (9) it is obvious that only the DM and u/d quark interactions 
break isospin effectively. Immediately, we draw the conclusion: if scalar interactions account 
for CoGeNT/DAMA experiments, the DM and first-family quark interactions must give the 
predominant contribution. Then the ISV /„ = Xfp transformed into quark level is simply 
given by 

a, - Bn - TBI ^ ^' ^'^^ 

and the ratio is about —0.77 for X = —0.7 (throughout this work, we shall use this ISV as 
a referred value, as well the referred DM mass Mdm = 8 GeV). Furthermore, the general 
effective DM-proton couplings with ISV can be organized in such a form 

It is factorized into DM-quark effective coupling and a model-independent factor casted in 
the bracket, which takes value —2.5 for I ^ —0.7. We typically require DM-proton SI cross 
section ap = 0.01 pb, in turn fp is determined to be around the value given above. By the 
way, it applies to both scalar and vector interactions in a notation given in Eq. (Bl). 

The above conclusion excludes the models with conventional Higgs mediator, e.g., in the 
one/two Higgs doublet SM and in the Minimal Supersymmetric Standard Model (MSSM) 
where the Higgs fields not only dominantly mediate interactions but also generate the SM 
fermion masses. In such models, we have 

rriui/d, Vu/d 

with Vu/d the vacuum expectation value (VEV) of Higgs field H^/d- It is independent on 
the mediator-quark couplings, consequently the second and third families give main contri- 
butions by virtue of larger form factor, thus the isospin is preserved. But if the Higgs field 
is a spectator Higgs field whose VEV is zero or very small, its Yukawa couplings with SM 
fermions are free parameters. And then the light quark contributions to SI crosse section 
can definitely exceed the contributions from heavy quarks. 

Now turn to the vector interactions. The quark bilinear matrix elements in the nucleon 
states are greatly simplified by virtue of the conservation of the vector current, to which 
the sea quarks and gluons do not contribute. As a consequence, the effective interactions 
between DM and nucleons are simply given by 

fp = bp = 2bu + bd, fn = bn = 2bd + bu. (13) 
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Notice that even in the case where these vector interactions are generated by other interac- 
tions, e.g., colored mediators discussed later, the above description still holds. So the ISV 
hn = 1-hp corresponds to the coupling ratio at the quark level 

^ = 1^ (14) 

for X = —0.7 it takes a value —9 : 8. Obviously, similar to the case for the scalar interactions, 
the ISV from vector interactions also must come from the interactions between mediator/DM 
and the first-family quarks. And the DM-proton effective coupling is 

3 

and fp = -1.25bd for X = -0.7. 

As given in the introduction, we should pay special attention on the inelastic DM 
(iDM) [9]. In the U{l)x model, the vector currents involving DM couple to the vector 
gauge boson non-diagonal, which in the mass basis are 

gabXal^XbZ'^, gab^id^^bZ'^, (16) 

where gab ~ Sab, and x is pseudo Dirac while (p is approximate complex scalar. In other 
words, for the ISV iDM, an f/(l)x niodel should receive special attention and we shall 
consider it in details later. 

Comments are in orders: (i) The scalar interaction realization of proper ISV depends 
on both the isospin- violations in the nucleons itself, i.e., i?|f^ ^ Buj^, and the DM-quark 
interactions, whereas the vector interactions only depends on the latter, (ii) The thermal 
DM with right relic density has annihilation rate {(Tanv) ^ 1 pb ~ lO^a^, given ap ~ 10"^ 
pb. On the other hand, due to crossing symmetry, the annihilation rate and scattering rate 
from the operators accounting for the direct detection experiments scale as follows 

—■ (I'J 

ap /i^ 



p 



It just gives the ratio at right order for LDM around 10 GeV. Interesting, this numerical 
coincidence involving three basic elements of DM, the mass, relic density as well as scattering 
rate. In this regard, it "justifies" the GoGeNT/DAMA inspired LDM models with ISV. Of 
course, the rough argument ignores the velocity suppression for the operators like XXQQ, but 
it is still sensible. Since in a complete model the pseudo-scalar type operators (x75X)(^75Q') 
usually do exist, and has comparable coefficients. And we will see it via concrete discussions 
later. 



B. Some Constraints 

For a light DM with mass around 8 GeV and having quite large DM-nucleon scattering 
rate, it suffers from a list of model independent constraints (denoted as Ci_4 in the following), 
coming from the cosmology, astrophysics and collider as follows 

Cii PAMELA The PAMELA measures the anti-proton spectrum from 1-100 GeV, and it 
has no deviation from background [17]. Since the LDM couples to quark substantially, 
it probably renders the low energy spectrum of the anti-proton excess due to crossing 
symmetry [18, 19]. This constraint holds when the LDM has rather large annihilation 
rate (> 0.1 pb) into quarks today, thus invalid for the DM annihilation with velocity 
suppressing. Moreover, such constraints usually can be avoided by choosing proper 
astrophysical parameters [18]. 

€2'- Sun Neutrino Because of the large DM-nucleon cross section, the Sun captures DM 
particles with large rate. These DM particles subsequently (cascade) annihilate into 
neutrinos, which are expected to be observed by the Super-Kamiokande. It excludes 
the thermal DM with predominantly annihilation modes into tt/uu/At as well as 
heavy quarks modes bb,cc [20]. We have to emphasize that such neutrino constraint 
is so strong that generically one has to sufficiently suppress the DM annihilation via 
these channels, especially directly to neutrinos. 

C-y- CMB The lack of distortions in the Cosmic Microwave Background (CMB) spectrum 
due to DM annihilations at redshifts z ~ 500 — 1000 may give constraint on the GeV 
scale DM [21]. To make DM being a thermal relic, its annihilation modes should be 
dominated by the /i or r modes (and corresponding 4 leptons), alternatively the DM 
annihilation is velocity dependent. 

C4: Colliders DM-nucleus recoil can be converted to the DM production at the hadronic 
colliders, and the lack of relevant signals provides another constraint [22, 23]. In the 
ISV scenario the Tevatron gives very strong constraint. Vector interactions (as weak 
as cTp ~ 0.001 pb) have been excluded definitely while the scalar interactions also 
close to the exclusion line, as given in the Fig. 1. It implies that the iDM with ISV 
must not be viable. To evade such Tevatron constraint from the effective operator 
analyses, we require that the mediator is light and its mass is much smaller than DM 
mass. In this way, the DM-nucleon cross section is enhanced by lighter mediator mass 
since o"„ oc l/^mediator whereas DM production at the collider scales as 1/s, so the 
constraint made in Ref. [22, 23] is evaded. In addition, Ref. [24] gives the LEP bound 



on the operators O-oueTe, and shows that this operator is not allowed to provide main 
annihilation channel of DM freeze-out. 

Let us summarize the constraints and point their possible implications (a recent more quan- 
titative similar consideration see Ref. [25]). Among them, Ci can be satisfied thus not very 
serious. But C2 is a universal constraint independent on the DM spin, and from it we may 
expect that DM annihilates into /i or light quarks since the final state e^ is excluded by C4. 
It is consistent with LDM in the CoGeNT/DAMA region where LDM likely annihilates into 
up and down quarks with cross section around Ipb, e.g., in the models given in Section IIIB 
and Section IIIC Then C3 may disfavor it. But take the astrophysical uncertainty into 
account, we just regard it as a referred constraint. Last but not the least, C4 gives the 
most powerful constraint in actual model building, it picks out the models whose DM- quark 
interaction either is a scalar type or via a light vector mediator exchange. In the case of 
light mediators as in Section III A, there is an elegant way to avoid various constraints by 
virtue of DM dominant invisible annihilation modes (to light states), and we shall discuss 
it in detail in the Section IV. 

III. LIGHT DARK MATTER MODELS WITH ISV 

In this Section, we systematically construct the ISV models. Let us explain our convention 
first. The s— channel Feymann diagrams have the vertices between DM-DM- mediators and 
between quark-quark-mediators, while the t— channel Feymann diagrams have the vertices 
betweem DM-quark-mediators. Thus, our s— channel and t— channel Feymann diagrams 
correspond to the t— channel and s/n— channel Feymann diagrams in the direct detection 
experiments, respectively. For simplicity, we shall not consider the DM models where the 
gauge bosons are DM candidates since they are generically complicated. Thus, the dark 
matter particles are either a scalar or a fermion. For s— channel Feymann diagrams, the 
mediators must be a scalar or a gauge boson since they can couple to two quarks. For 
t— channel Feymann diagrams, the mediators must be a color triplet fermion if DM is a 
scalar. And the mediator must be a color triplet scalar or vector boson if DM is a fermion. 
Because the models with color triplet vector boson are complicated in general, we will not 
study them in this paper. Therefore, we shall consider three kinds of LDM models with 
ISV, where the mediators are dominantly the U{l)x gauge boson Z', a spectator Higgs 
doublet, and color triplets respectively. For the first and second kinds of models, mediators 
propagate in the s— channel, thereby ISV is independent on the DM-mediator couplings, 
which simplifies the discussion. However, the color triplets propagate in t— channel, and 
their analyses are a little bit involved. 
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I CoGeNT mod 
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FIG. 1: Top: favored regions and exclusion contours in the (Mdm, o"n) plane, ISV frj fp = —0.7, 5 = 
0; Bottom (for the annual modulation): ISV fn/ fp = —0.7, 5 = 15 keV. Collider exclusion for some 
operators are also imposed. The data are from Ref. [10] and Ref. [23]. 

A. Chiral U{l)x Vector Boson 

A simple way to produce ISV is to introduce an exotic U{l)x gauge boson as the mediator. 
From the discussions in the last Section, we only consider the light Z' . In practice, some of 
the following discussions that merely involve in the direct detection apply to any Z' boson 
with mass much larger than the transfer momentum. First, we briefly prove that the Z boson 
in the SM can only generate ISV with \fn/ fp\ ^ 1- Explicitly, the Lagrangian between the 
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Z boson and SM quarks 
1 



(^NC ^ 



COS Q,, 






1 



2 3 

1 1 



sin^ 9yj]uL + ur'j'' 



sin^ 9^ ur 



sin^ 6u,]dL + dnY 



2 3 
Utilizing the formula Eq. (13), one can easily get 

for, 



sin^ 6,,, 1 d 



R 



(18) 



(l-4sin=^ft^) ^ -0.08 < 1. 



(19) 



We emphasize that this result depends only on the SM structure. 

Now consider Z' from an exotic U{l)x- We can gain some insight into the structure 
of the quark charges under U{l)x, by investigating the ISV origin. And the discussions 
have nothing to do with concrete information in dark sector. We start from the Z' and SM 
fermion current couplings 



CnC = -9X Y. f^^' [iQf.)^^PL + [QfnhPR] I)K^ 



(20) 



«j 



where gx is the gauge coupling of U{l)x and Q/^/^ are the charge matrices of the left- and 
right-handed fermions in family space. To not induce the tree-level flavor changing neutral 
current (FCNC), we reasonably assume they are diagonal matrices (for Abelian gauge group 
it must be ture). In the quark sector, transforming quarks into the mass eigenstates via 
Ql/r -^ ^ql/r'^l/R (^^^ same letters are used to label gauge basis and mass basis), we have 






NC 



q=u,d;i,j 



q,i 



(21) 



In the second line, we explicitly show the first family, while dots collects axial vector current 
and other families. 

In this paper, we consider that the charge matrices are family universe {Qq^/^i = Qqi^mi 
then the unitary of Vg^ ,^ ensures that the SM fermion mixings will not induce flavor violation 
in Eq. (21). And the first-family quark charges give rise to proper ISV: 



K 
hd 



^ul ' 'V«H ^qL ' ^ 

Vdi ' ^(Ir ^qL 



^UR 



(22) 



In the second equation QuL=QdL= Qq^ is taken since ul and d^ are in the same SU{2)l 
multiplet. Clearly, the ISV effect is ascribed to the different right-handed up and down 
quark charges. So quarks must belong to the chiral representation under U{l)x, and we are 
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forced to work in the two Higgs doublet model (2HDM) where Hu couples to up-type quarks, 
while Hfi only couples to the down-type quarks and charged leptons. By the way, it is easy 
to check that such U{l)x can not come from Eq gauge symmetry if we want X ^ —0.7. 

Through simple arguments, some useful conclusions are draw. The U{l)x gauge invari- 
ance of Yukawa interactions give 

Q,, - Qun + Qh^ = 0, g,, - Qdn + Qh, = Q. (23) 

If we do not introduce additional colored vector-like fermions, the 5'[/(3)^?7(l)x anomaly 
cancellation condition gives 

3(q,.-%^-^)=0. (24) 

From these two equations, we get Qh^ = —Qh^, and then we just need to introduce one 
Higgs doublet for the non-supersymmetric models. Combining them and Eq. (13), we get 
the charge conditions to obtain proper ISV 

_7-5X n - ^~'^^ n n - ^^^ n 

And X=-0.7 leads to a somewhat peculiar solution: Qup, = 35/34(5 h„, Qd^ = — 33/34(5 _?/„ 
and Qqj^ = 1/34(5 h„ [44]. However, if Eq. (24) is relaxed by introducing new colored vector- 
like fermions, the elegant solution may be found. But we just assume Eq. (24) is held in 
this paper. 

Finally, the large DM-nucleon recoil cross section is easily achieved with light mediator. 
Assume the Dirac or complex scalar DM carries U{l)x charge Qbm, then in light of Eq. (15) 
and Eq. (25), we have 

r _ QxQ'dmQhu 2 

Jp ~ ' 



M|, 1 - X 



where X = —0.7 is took. So typically we require a very weakly coupled (at least to quarks) 

U{l)x. 

Furthermore, the DM particles can be a scalar or a fermion. In particular, the U{l)x 
gauge symmetry can be broken down to a Z2 discrete symmetry, which stabilizes the DM 
particles. Thus, we do not need to introduce extra discrete symmetry for DM particle to be 
stable. 
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B. (Approximate) Spectator Higgs Doublet 

Higgs doublet is another typical mediator. But as argued in Section II A, if the Higgs 
doublet mediates DM-quark interactions meanwhile accounts for the SM fermion masses, 
the allowed ISV will be ignorable. So we have to introduce the Higgs doublet that is a 
(approximate) spectator to the electroweak symmetry breaking, namely it is only in charge 
of mediating the DM-SM fermion interactions. The origin of ISV arises from the Yukawa 
interactions involving in quark sector 

-Ly D yyu,iqL(^HluR + yd,iqLHidRJ + (yu,2qL^HluR + yd,2qLH2dRJ + h.c. . (27) 

Here Hi is the conventional SM Higgs fields whose neutral component acquires a VEV 
{Hi) = V around 174 GeV. While H2 = (h+, ^^^^^) , which carries the same SM 
quantum numbers as Hi, is the spectator Higgs with ignorable VEV. Thus, yu/d,2 are free 
parameters, which we shall address this issue in details later. Quite interestingly, recently 
such a spectator Higgs doublet is also inspired to explain the CDF W+2jets anomaly [35, 36]. 
Now we shall consider the spectator Higgs and dark sector interactions. To keep the 
discussion as general as possible, we do not specify to any concrete model. The CP-even 
particle if ° mediates DM-nucleon SI scattering, and as usual, H^ probes the dark sector via 
the operators 

al^pi/o, x(« - /37^)xi7°, (28) 

a^l^pA", xiaA-f3Al')xA', (29) 

where the CP-odd Higgs field A^ does not generate SI cross section but opens another DM 
annihilation channels. In our scenario, we assume that H^ barely mixes with the CP-even 
Higgs component of Hi, thus, such mixing factor can be ignored. The DM particles can be 
stable by introducing the Z2 disccrete symmetry. If the DM particle is a real or complex 
scalar through Higgs port [38], we can realize the first term in Eq. (28) easily via the following 
renormalizable term 

X^\(f)\'^HiH2 + h.c.^a = V2X^v. (30) 

If the dark matter is a fermion, to have the renormalizable interactions between the DM and 
H^, one may consider a dark sector containing term XSxx^ aiid the singlet scalar S further 
mixes with H^ after the electroweak symmetry breaking. A more interesting possibility arises 
in the supersymmetric models, e.g., in the next to the MSSM where the light neutralino xi 
is the DM [39]. 
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The ISV is determined by the flavor structure of H2 Yukawa couplings with quarks. 
Transforming quarks into the mass eigenstates, we find that Eq. (27) turns to be 

Cy D /= uYuuPru - H'dYduPRU 

+ y=—dYddPBd + H+uYudPRd + h.c, (31) 

where the effective Yukawa coupling matrices are defined by Yqqi = V^^yq'^2Vq'^R. Generally 
speaking, F— matrices are not diagonal, consequently the spectator neutral Higgs mediates 
tree-level FCNC. But yq^2 and Vg^ are free, thus, the flavor problem may be solved by choosing 
them properly. Integrating out H^, the operator coefficients for the SI cross sections are 



a, 



q 



V2(2Mdm) ^m '' '' ' V2'. 



m. 



"HO 

And in light of Eq. (10), the ISV X ^ —0.7 is obtained provided that the Yukawa coupling 
ratio (Yuu) 11 /(Ydd) 11 ~ —0.77. Of course, we also require (^^5)22,33 sufficiently small so that 
the DM interact dominantly with the first-family quarks mediated by the scalar interaction. 
Now we further assume the spectator Higgs is the source of the W-|-2jets anomaly via 
the process pp — )■ H^ — )■ W^H^/A^ — )■ (i^(i~u +33 [35]. As a consequence the measured 
invariant mass of two jets fixes the charged Higgs mass m^^o ^ 150 GeV. For the benchmark 
scenario given in Ref. [35]: {Yuu)ii — 0.06 and mH± ~ 250 GeV, then to produce Up ~ 0.01 
pb, we obtain 

a (F„„),, ~ 1.6 GeV, a {Ym)^^ ~ -2.0 GeV, 

a (KJn ~ 0.10, a (Fdd)ii ~ -0.13 . (33) 

For a real scalar or Majorana DM we only need half of these values. But H^ will decay to 
a pair of DM, with branch decay width ratios over the one of H^ — )■ uu as follows 

r(g° ^ 2DM) (l^Jr/ ( a{yuu)ii V 1 ( a{Yuu)ii V . . 

T{H^ ^uu) ^ 3 \ Mh, ) ' ?>\ iY^u)\, ) ' ^ ' 

In the above estimation we have used the tree-level decay for H^ decays to quarks, nev- 
ertheless the QCD corrections are important. For the (real) scalar DM the two chan- 
nels are comparable, thus, the scenario is viable. Whereas for the fermionic DM H^ pre- 
dominantly invisible decays to DM. The situation can be improved by considering process 
pp — )■ H*^ — )■ W^H^ — )■ £"*"£" z/ + jj [35], further inverting H^ and H^ mass hierarchy. Then 
H^ — )■ H^^+DM+DM decay width is suppressed by additional phase factor factor l/(27r)'^. 
Moderately increasing {Yuu)ii and/or lowering rriHO so as to increase the production rate of 
if^, the fermionic DM may also be consistent with W + jj. 
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Finally we investigate the condition that such (quasi) spectator Higgs doublet appears. 
Equivalently speaking, when V2 is small enough, the Yukawa couplings yu/d,2 become free 
parameters. Roughly it requires (Kn)ii^'2 < "^n, (Xid)ii^2 < "i^, where {Yuu)ii ~ 0.06 
are taken to explain W+2jets. Thus, we have V2 < 0.05 GeV. Now let us consider the 
renomalizable scalar potential of two Higgs doublet model 

V =^A\H^\' + ^ll\H2? + ^iUh\h2 + hIh{) + ^|i7i|^ + ^iH^t 

+ A3|i/ir|i/2p + A4|i/|i/2p + ^{HlH.Y + ^{HlH,)\ (35) 

where all parameters are assumed to be real for simplicity. The small V2 can be achieved by 
a small Higgs mixing term /i^2 ^ 0(10) GeV^, which can be seen by the tadpole equations 

dV 

-—— = HI2V + A2W2 + hv'^V2 + \4v'^V2 + Asf ^t;2 = 0, (36) 

oho 

from which it gives V2 /i?2^/((-^3 + A4 + As)^^ + /x^) ~ O(10)MeV with A3,4,5 ~ 0(0. 1). 

In conclusion, it is not difficult to get an approximate spectator Higgs doublet to account 
for ISV as well the W+jj anomaly, given a scalar DM. 



C. Color Triplets 

The color triplet mediators are distinguished from the Z' and Higgs mediators since 
they mediate DM-quark interaction in the t— channel. In addition, they usually generate 
hybrid interactions, namely generating scalar and vector interactions simultaneously. It is 
more convenient to study such models from complete Lagrangian directly. In this paper, we 
concentrate on the fermionic DM particle and scalar color triplet mediators. The discussions 
for the scalar DM particle and fermionic color triplet mediators are similar, and will not be 
studied here. Moreover, to stabize the DM particles, we consider the discrete Z2 symmetry. 
In particular, the DM particle and the color triplet mediators should be Z2 odd. 

We consider the most general mediators structure, each quark type Ul/r and di/R is 
accompanied by a color triplet mediator Ql/r, just like the corresponding squarks in the 
MSSM. Since we are aiming at ISV, only the first family will be considered. Then the model 
preserving SU{3)c x f/(l)em gauge symmetry takes a general form, adopting notation easily 
matches to the MSSM: 

^ = - ^xXX - J2 "^L/3^°4 ~ Yl ^d^p^^d^l3 

a,l3=L/R 
a 

- E [^^-^^(l + ^5)951 + A;^x(1 - l5)qql + h.c] . (37) 

q=u,d 
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Since qi and q^ carry identical 5'f/(3)c x U{1)em quantum numbers, we allow triplets 
left-right (L-R) mixing as well as quark-mediators coupling with opposite chirality (L-R 
coupling), namely Ag^ ^ and A'^ 7^ 0. They are sources of chiral symmetry breaking. The 
triplet mass eigenstates are denoted by gi 2, related to gauge eigenstates by q^ = Ylii ^qMi'- 



Fq^l = COS 6*5, 


Fq^2 = sin 9q, Fq^l = 


-sin^g, Fq^2 = 


= COS^g. 


(38) 


an Uq ^Xq -I 


/l + x2 < 0, Xq = (m|^ 


(39) 



The corresponding mass eigenvalues are given by rriq-^^ ^ . Then the interactions can be rewrit- 
ten in a from 

£ D - X ("f + I3?i5) qqi + h.c, 

Z. (^^ + KJ ^..'' A' = E i^^- - KJ F'^^i- (40) 






In light of Eq. (A9) and Eq. (A12), integration out $ leads to the effective operators 
involving SI cross section in the form of Eq. (2) and Eq. (3). And the operator coefficients 



m 






bu = — 

Also, the expressions for a^ and 6^ are obtained by replacing u with d in the above equations. 
There are two interesting limits. One is the chiral limit a^jd — )■ 0, arises when both L-R 
mixing and L-R coupling are ignorable, then the scalar interactions vanish [45], leaving the 
pure vector interactions. In other words, interactions mediated by triplets in the chiral limit, 
the iPM scenario can be realized even in the absence of a vector gauge boson . This novel 
phenomena should receive special attention, but such models can not give negative X since 
hu and h^ take the same sign. In contrast, in the Majorana limit where the DM is a Majorana 
fermion (or real scalar), the vector interactions vanish. In conclusion, colore triplets must 
generate ISV via the scalar interaction, and chiral limit must be avoided. In general it is 
expected that we both the scalar and vector interactions, as well as /« = «„ -|- 6„. 

L-R couplings usually are small, since the SM quarks and corresponding mediators are 
chiral and then are distinguished. In the MSSM, the lightest supersymmetric particle (LSP) 
is usual neutralino. The Higgsino component in the LSP neutralino has such coupling from 
superpotential, which nevertheless is tied to the chiral symmetry breaking in SM and then is 
tiny. Of course, one can relax this constraint by recurring to neutralino from spectator Higgs 
doublet, but the following discussion will not much different to the previous subsection, so 
we do not take it into account and only consider the L-R mixing as unique chiral symmetry 
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breaking source. So we have A^^ = A' = 0, and 



9l 'to 

a, = -^ X (!!!3fl!!k,i^20,] ^ -^F„ (42) 

with q = u,d. Clearly, if the two states are highly degenerate or L-R mixing is small, the 
effective coupling will be suppressed. So we consider the general case with Fg ~ 0{1), as a 
result, according to Eq. (11) with X = —0.7 

for complex scalar double this estimated value. 

An application of this result is the supersymmetric model with light t/(l)x— gaugino LSP, 
and the first-family squarks mediate DM-quark interactions (the second and third family 
contributions can be suppressed due to heavier masses or smaller U{l)x charges). Relevant 
terms are 

/2 ~ 

CD- —Qq^gxBx (1 - 75) qql 

/2 ~ /2 ~ 

- -YQ^r9xBx (1 + 75) uv}r - -^QdRgxEx (1 + 75) c?4? + ^-c-, (44) 

where Q denotes U{l)x charge. The couplings defined in Eq. (37) are inherited from gauge 
couplings 






^dL -^UL - -rrdxQqL^ ^dR - -7r9xQdR, A„^ - -^gxQun- (45) 



But actually a^ are only semi-quantized, by virtue of extra parameters Fq/m'~_^. And now 
we have 



ad QdR \mlFd 



-0.77. (46) 



From Eq. (39), it is seen that Qua/QdR <0 is still a necessary condition to get X < 0, 
implying a chiral U{\)x- 

In practice, the MSSM with U{1)y bino dominated LSP just fits this scenario. However, 
it fails to explain CoGeNT/DAMA quantitatively. In light of Eq. (43) and Eq. (45), even 
we have Fd ~ 1, to make squark masses lie above collider lower bound (^500 GeV), we 
still need large quark gauge coupling gxQq ~ C'(l) (not excluded by collider for sufficiently 
heavy Z'). Whereas the (s)quark f/(l)y gauge couplings are too small. We have to stress 
that, exotic U{l)x is still interesting enough, since it realizes ISV in an economic way, asides 
from introducing U{l)x vector multiplet. Moreover, it does not trigger flavor problem. 
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Finally comment on another interesting aspect of the existence of color triplets at low 
energy based on SM. In the model given by Eq. (37), the Yukawa couplings Xq/u/d generically 
are complex c-number, so they introduce new physical CP phases. Recall that the SM fails to 
give sufficient baryon asymmetry, then these new color triplet mediators (If they are scalars, 
no anomalies will be introduced.) not only account for proper ISV, but also potentially 
provides adequate enough baryon asymmetry via heavy color triplet decays. Amazingly, for 
a Dirac dark matter, it is possible to realize asymmetric DM by the triplet common decay 
so as to resolve the coincidence puzzle ^dm^^ : VLi,h'^ ^5:1. By the way, non-annihilating 
asymmetric Dirac DM automatically satisfies astrophysical constraints. We believe it is a 
minimal unified framework to realize ISV, asymmetric DM as well as baryogenesis, and we 
leave it for future work [37]. 

D. Dual Mediators 

In the previous discussions we concentrate on ISV arises due to a single-type mediator, and 
no conventional mediators succeed in giving proper ISV. However, as proved in Section III A, 
the SM Z— boson actually only mediates DM-neutron interaction. So combining Z boson 
with Higgs field (without ISV effect) in principle can produce proper ISV. In fact, a complex 
sneutrino DM vi [30, 31] just fits this scenario. Unfortunately, Z mediator plays no role in 
the ISV scenario. To satisfy the constraint from the Z— boson invisible decay (into a pair of 
DMs) width [33], we obtain the model independent constraint on the vertex gz^^Z^'-Uid^i'l, 
whose effective coupling must be smaller than 0.023, i.e., gz^ < 0.023. Then we have 
/„ ~ — M^'^g gzj^j^/ 4: cos 9w ~ —0.5 x 10~^GeV~^, far to enough in the ISV scenario. 

Another mediator frequently appears in the literatures is an exotic X^ boson of U{l)x 
under which the SM fields are neutral, whereas U{l)x and f/(l)y have kinetic mixing term 
so that Xfj^ can mediate DM-quark interaction. As noticed in the Ref . [7] , this kind of X^ 
only mediates DM-proton interactions. Let us show this point explicitly. We consider the 
gauge kinetic sector oi U{l)x x f/(l)y 

t-gauge = ~'^^Y ^Yfiu ~ 'l^X ^Xfii/ + ^-^y rXfiu ■ (47j 

We are interested in the small mixing limit ^ ^ 1, since it will not significantly affect the 
Z invisible decay width even with U{l)x charged LDM. Moreover, we consider the collider 
favored light X^ scenario therefore Mx < Mdm <^ Mz. Then, after the electroweak symme- 
try breaking, working in the mass eigenstate basis {Zf^,A^,X^), we obtain the interactions 
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between the X^ and the SM current at the leading order of 6 [34] 



Jem =9 sin 9n 



-uY'u + ( — ) d-f^'d + (-1) e-f^e 



(48) 



Apparently, the kinetic mixing just induces the coupling between the X^ and electromagnetic 
current J^^. So X^^ behaves like the photon that also only mediates DM-proton interaction. 
We have shown that X can violate isospin. Next we consider the DM and light gauge 
boson gauge interactions, including the following DM-X^ effective couplings 

-Cdm 3 ^DM^M0*^^0, gouX^xYx- (49) 

So the DM-proton interaction effective operator coefficients are given by 

fg sin2e^\ /^dm^\ /^„x 

The X^ alone is not adequate to produce right X, so we introduce the conventional Higgs 
as another mediator. The relevant discussion has been represented in Section IIIB, here 
we just identify H'^ as the SM Ifiggs h, then the DM-quark effective couplings are given by 
Eq. (32). In this case in the Eq. (7) ag/niq oc 1/v is family universal, and we get the isospin 
conserved contribution 

^n Q 1 ( f{n) n ^ („) 

"^""""v^(2M,m)^<V"^ +'"27^-- 

_5/10GeV\ /lOOGeVV ^ ,,_2 



= 0.15 X 10^^ — GeV"'. (51) 

V Mdm ) \ nih ) ^ ' 

We only include the s quark and heavy quarks contribution that gives a value ~ 0.31 in 
the bracket. Moreover, we assume that the A^ji/pl^p term gives the coupling constant a 
as well as A^ = 1. Perturbitivity and naturalness (large A,^ contribute sub-TeV mass to the 
DM) do not favor large A^. But the two Higgs DM model with large tan/3 enhancement 
can solve this numerical problem. For example, consider a simplified case where couples 
to Higgs fields mainly via a term \^\(l)\'^HdHu, then in the Eq. (51) via H^ mediation (we 
assume the heavier CP-even Higgs is dominated by if^) 

-^^^?^ = v^A,tan/3. (52) 

V Vd 

and now it is required that A0tan/3/(mj|^o/lOO GeV)^ ~ 1. In conclusion, the total DM- 
nucleon coupling fp = ap + bp,fn = an^ap, and the proper ISV is produced as follows 
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which is readily to be reahzed. 

Analysis on fermionic DM can be employed similarly, with the replacement a/2M£)M — ^ ot- 
In general, the numerical problem is exacerbated, and maybe somewhat light vnh (even with 
large tan/3) is required to enhance the cross section, just as required in the NMSSM with 
neutralino LSP [39]. 

IV. SNEUTRINO IDM WITH LIGHT Z' MEDIATOR 

Interestingly, the light iDM models with ISV can explain the CoGeNT and DAMA annual 
modulation consistently, as well as reconcile them with all other results [10]. In particular, 
we can evade the other astrophysical constraints naturally in quite a few models, for example, 
the spectator Higgs mediator models. To be concrete, we shall propose sneutrino iDM model 
with a very light and very weakly coupled Z' mediator. 

A. Motivation for iDM 

In the first place, in light of the general analyses made in Section II B, the Z' should be 
lighter than the DM. We even require Mz' < 1 GeV to escape from various astrophysical 
exclusion elegantly. It is much lighter than the consideration like in Ref. [27], nevertheless 
natural since the Z' mass scale closes to the LDM mass scale [46]. The light Z' mediator 
allows very weakly coupling gf^ ^ 1? in this way the Tevatron [22, 23] and LEPII [40] con- 
straints are satisfied. Moreover, kinetic block mechanism makes the DM main annihilation 
products to be soft 2e or/and 2/i today. Then the constraints from Sun neutrino, CMB as 
well as PAMELA can be readily satisfied. 

We explicitly show how this mechanism works for a scalar iDM. At the early Universe, the 
DM has two kinds of comparable annihilation channels: DM+DM— )■ qq via Z' mediation, 
as well as the invisible annihilation modes DM+DM—)- XX, with X the sub-GeV hidden 
states such as Z' (we take it as an example). Then the total annihilation rate is 

where Qf is understood to be effective charge in the mass eigenstates. Terms in the bracket 
collect DM annihilation into hidden states. Provided that the U{l)x charges for the SM 
fermions and DM particle are comparable, today the second term overwhelmingly dominates 
(since v ~ 10"'^), subsequently the Z' cascade decays into e//i/7 have no dangerous modes, 
thus the strong tension with astrophysics is resolved. 
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In contrast, the Dirac iDM annihilation modes DM+DM— > // have no v'^ suppressing, 
so it is quite unnatural to suppress dangerous modes to desired level. In fact, direct an- 
nihilations into the SM fermions (assuming family universal Z') are almost fixed by cross 
symmetry. According to Eq. (26), to get the adequately large DM-proton (inelastic) scat- 
tering rate ap ~ (9(0.1) pb, we need q'xQh^Q-dm ~ lO"'^. On the other hand, from Eq. (54) 
we know that the direct DM annihilations into the SM fermions have a rate about (9(10~^) 
pb. While the Dirac DM annihilations via s— wave have even larger rate (9(10^^) pb. This 
result is derived from cross symmetry, thus, the Dirac DM annihilation rate can be reduced 
by no way, except by lowering Z' mass even orders lighter. Or one can consider only the u 
and d quark charges under U{\)x to avoid Sun neutrino constraint, but it leads to a weird 
model. In conclusion, scalar iDM with sub-GeV Z' mediator moreover dominantly decays 
to sub-GeV light hidden states is a preferred model for the GoGeNT annual modulation. 

B. Light Sneutrino iDM from Low Scale Seesaw Mechanism 

In the MSSM extended with low-scale seesaw mechanism, (light) sneutrino is a natu- 
ral scalar iDM candidate (We also present a quite simple non-supersymmetric model in 
Appendix C). To realize ISV we further consider a light f/(l)x— extension, and the right- 
handed neutrino (RHN) is charged under U{l)x- Then it is natural to attribute the origin 
of Majorana mass scale to the U{l)x breaking scale. The relevant parts in the minimal 
model are 



W Dy^^L,N,H^ + -^SN^ + fiH^H,, (55) 

-Csoft 3 (m^JiV.p + m||5p) + Ao Lf^L^N^H^ + ^SN^ + h.c\ , (56) 

For simplicity, all parameters are assumed to be real, and soft mass square terms for the 
MSSM singlets S and A'^j (we simply consider one family only) are around GeV scale or 
even smaller. To reduce the parameters, we assume a common trilinear term Aq. Because 
it controls the mass splitting of A^j, we further assume it to be very small on our purpose. 
Such soft parameter pattern is natural in the gauge mediated supersymmetry breaking 
(GMSB) scenario. S carries U{l)x charge Qs = —Qn/'^, and it develops a VEV with 
(5) = Vg ~ O(IOO) GeV. Thus, it breaks the exotic U{l)x gauge symmetry at low scale and 
gives the Majorana mass terms for the RHN M]si = Xvg. Since we want to naturally have a 
sneutrino state with mass around 8 GeV to be the LSP, we expect M^ ~ 0(10) GeV, and 
then the neutrino Yukawa couplings y'^ ~ (9(10~^) are irrelevant to our discussion. 

This model is able to realize the above scenario readily. First, let us consider the f/(l)x 
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spontaneously breaking. The scalar potential of the S is given by 



Vs =Vf + Vd + Vsoft, 



VF=\XiSNi + y'iLjHu\', 



Vn=^(QHjH^\' + QHjH,\' + Qs\S\' + Y,Qf\fn ' (57) 

and Vsoft = —^soft is given by Eq. (56). So the crucial part of the scalar potential of S can 
be casted in a form of simple (j)'^ field theory 

Vs^-^^l\S\' + l\S\\ (58) 

where 

^^l = [qIQh^Qs cos 2/3) v^ -ml, k = 2g\Q% (59) 



where Qh^ = —Qhu is used. Interestingly, a simple way to trigger f/(l)x breaking is available 
even with positive m|: the first term in the /i| is negative provided 



Qh^Qs < 0, (60) 

then /x| in principle can be as (positive) small as will. It leads to 

Vs = y-f^s ~ ^/cos 2 (3\QhJQs\ X V, (61) 

at the weak scale. In the above estimation we have set ml — )■ 0. In other words, in our 
model D— term is adequate to break the U{l)x gauge symmetry without resorting to any 
exotic requirements. 

Next we examine the mass spectrum. The Z' mass is determined by Vg- We have to 
stress that without S"s VEV one linear combination of Z' and Z.. is still massless after the 
electroweak symmetry breaking. The point is that Qh^ = ~QHa and then the Higgs mixing 
term B^HuHa is allowed. Now, ignoring the small Z — Z' mixing effect, at the leading order 
we have 



Mz' ^V2gx\Qs\vs ^ V2 1 g^J I g5 1 cos 2/3 x gxv . (62) 

Thus, Mz> can naturally at the GeV order by choosing gx ~ 10~^. On the other hand, at 
the leading order, the CP-even Higgs field from the singlet S (denoted as hg) has degenerate 
mass with the Z'. However, in this minimal model, LSP tends to be singlino 5", which has a 
light (degenerate with Mz>) Dirac mass term with the U{l)x gaugino A. So the heavier the 
Majorana gaugino mass term of A, the lighter the Majorana mass fermion S. To solve this 
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problem, we introduce extra singlets 5" (they are maybe required to cancel the f/(l)x gauge 
anomalies), with terms MiSSi or S'^Si. As a result no light (lighter than A^) i?— parity odd 
particles are left in the particle spectrum. 

In this model, DM has many annihilation channels, and we simply assume the dominant 
mode is DM+DM— t- 2hs via the contact interactions from |Ap|A^||S'p. Right DM relic 
density requires A ~ 0.1, which is consistent with Vg ^ O (100) GeV and Mjv = Xvg ~ 10 
GeV. 

Finally, the sneutrino LSP is dominated by the RHN sparticle N. Taking into account 
that the mass splitting between CP-even and CP-odd states due to non-zero Aq, we get an 
iDM with mass and splitting 

MlM = m%^ = >^'vl + m%, 5^^Ao. (63) 

So to get S ~ 10"^ GeV, Aq is expected to be around this scale. It is maybe upset by virtue 
of fine-tuning, but a version equipped with inverse seesaw mechanism [29] can resolve this 
tuning problem [30]. This model contains an extra singlet N, 

W D VnLHuN + niNNN + Mj^N^/2. (64) 

Then the sneutrino mass splitting can be related to the origin of small neutrino mass rather 
than merely the soft mass terms [30]: 

mTv ruN B 

r^ m„ l-rriy , (65) 

rriD m/5 m^ 

where we have assumed the Dirac mass term rriD = unVu is much lighter than the second 
Dirac mass term mj^: tti/) <^ m^^ ^ m^ . And B is the bilinear soft mass for A". Even set 
i? = 0, in principle mN/mo ~ 10^ just give the right order for 6. 

V. DISCUSSION AND CONCLUSION 

Light dark matter models with the ISV fn/fp ~ —0.70 and large DM-nucleon spin- 
independent cross section cr„ ~ (9(0.01) pb may provide a way to understand the confusing 
direct detection experimental results. Combining with stringent astrophysical and collider 
constraints, we can further deduce the DM properties. In this work, we investigated the 
possible origin of ISV. General anatomy on the effective operators generating SI cross section 
is made, and we found that ISV essentially arises from the DM and first-family quarks 
couplings. To further explore their UV origin, we propose three kinds of models with the 
following mediator structure: 
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Z' from f/(l)x The U(V)x must be chiral, and the hght Z' is strongly favored. 

Spectator Higgs Conventional Higgs doublet mediates interactions preserving isospin. So 
we have to introduce (approximate) spectator Higgs doublet whose couplings to the 
SM quarks are free parameters. Such a Higgs doublet can be used to explain Tevatron 
CDF W + jj anomaly. 

Color triplets Combining the squarks in the MSSM with exotic f/(l)x (quarks strongly 
charged under it), we found that the light Bx LSP can generate proper ISV via the 
first-family squark mediation. This scenario is economic furthermore suffers no flavor 
problem. 

Exotic Z' plus Higgs For a SM-neutral \J{\)x having kinetic mixing with ?7(l)y, its light 
gauge boson X^ only mediates DM-proton interaction. Combining it with a conven- 
tional Higgs mediator, we can obtain the right ISV. 

Inspired by the CoGeNT annual modulation, we propose a light sneutrino as the iDM. It 
is charged under a f/(l)x, and the light Z' [Mz' <1 GeV) leads to a proper ISV. Moreover, 
this model naturally reconciles the ISV scenario with the stringent constraints from the 
astrophysical (such as Sun neutrino) and collider (Tevertron mono-jet search) experiments. 
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Appendix A: Operators after Integrating Out Mediators 

The models on the DM-SM fermions interactions can be classified based on the propaga- 
tors mediating DM and SM particle interactions. In this appendix, we borrow some results 
from the Ref [13]. For the scalar DM, it interacts with SM fermions by exchanging a Z' 
boson, a (real) Higgs doublet h and the colored fermion Q, the corresponding Lagrangian is 
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given by 

£ = -If^^F;, + im|, Z'^'Z; + a0t9^0Z'^ + qY{a - Pl')q Z'^, (Al) 

£ = 1 (9/i)2 - im^ /i2 - acP^^h - q{a - f3-f^)qh, (A2) 

C = Qitd-mQ)Q-q{a- /37^) Q0t _ z^^.. (A3) 

Integrating out the heavy propagators via equation of motion, we obtain the effective oper- 
ators generating SI cross section (other operators belongs to Eq. (4), we do not list here) 

^eff 3 -^(0t9,0), (A4) 

Ceff D -^(l)^<pqq, , (A5) 

£e// D — (|«p - 1/3^) gg0V + ^ (l«r + l/3r) gT'^g^^ V (A6) 

mQ rriQ 

For the real scalar DM, vector interactions disappear. Interactions of fermionic DM can be 
described analogous to scalar DM 

C = -\r^^r^^ + \m\, Z'^Z\, + X7^(a - h')xZ' , + qi^ia - h')qZ',, (A7) 

C = ^{dhy - ^mlh' - x{a - h')xh - q{a - h')qh, (A8) 

C = \d^\'^ - m||$p - x{a - /37^)g$ - h.c. , (A9) 

where $ denotes the scalar color triplet mediators. And the corresponding effective operators 
are 

Ceff D ^aa xYx qi/^q, (AlO) 

aa _ _ 
'"'h 
Ceff ^ ^ [(!«!' - W XXqq + (l«P + I/SP) Xl.Xqrq] ■ (A12) 

When the DM is a Majorana fermion or a real scalar, the vector interaction vanishes. 



Appendix B: Scattering and Annihilating 

In this appendix we briefly introduce the formula involving direct detections and give the 
relevant annihilation rates. Ignoring small ISV from the form factor of proton and neutron, 
the DM-nucleus SI scattering cross section at zero momentum transfer (not the actual cross 
section) can be written in a form [41] 

a,= ^-^[Zf,+ {A-Z)U]\ (Bl) 

n 
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where A is the atomic mass of the nucleus while Z is its atomic number. The reduced mass 
fiN = Mdmit^'n/^Mdm + rriiy), and (5c = 4 for the self-conjugate particle like Majorana and 
real scalar DM, otherwise 5c = 1- In the reference DM-proton scattering cross section is 
used frequently 



a, 



5cp^. 



2 
P f2 



TT 



fi. (B2) 



where fip the DM-proton reduced mass. And fp from scalar interaction in Eq. (2) and vector 
interaction in Eq. (3) are respectively given by 

Fermionic DM : /^ = ^ E^a,; fp = 2a„ + a^, 

Scalar DM: f^ = y^BP-^; /^ = 2a, + a^. (B3) 

The DM can annihilate into the SM particles, which determines the final relic density after 
the annihilation freeze-out. And the thermal average annihilation cross section, expanded 
with relative velocity Vrei (subscript will be omitted) takes a form of 

{av)F.o. = a + h{v^) = (a + 36/x/) (B4) 

where xj = Mdm/Tj = 3/(t'^) with Tf the DM decoupling temperature, for the weakly 
interactive massive particle (WIMP) typically having xj ~ 20 — 30. The DM relic density 
can be formulated as 

^,2 1.07 X lO^x.GeV-i , , 

nh^ ^ —- — — --^ , B5 

Mpiy/g^{av}F.o. 

with g^ as the effective relativistic degree of freedom when the DM decouples. The actual 
effective annihilation rate, which is used to determine relic density and calculate signals from 
DM annihilation today, is given by 

{crv)F.o = Tbm X (a -h Sb/xf) , (B6) 

wher Tdm = 1/2 for the complex DM while Tdm = 4 for the self-conjugate DM [42]. 

The a and b can be extracted out from partial wave expansion of the cross section times 
the relative velocity av = a + bv"^ (regarded as a rough thermal averaged cross section). 
Here the relative velocity is w = 2^1 — AM'^y^/ s with s the Mandelstam variable. For the 
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fermionic DM, av from operators involved given by [43] 



afXXff ■ 
^Xl'xh'f : 

^xi'xlf . 

Gspj 



3 „,2 



xxfi^f ■ 



V2 
bfXl,^xh''f ■ 

^xi'i.xh'i'f ■■ 

"^Xl'l.xHf : 



16^ X 2«/^dm/3> 
X GgpjM^y^l3f, 



IGtt 



47r 



X 



26^M^m/3/ (2 + zf) 



^ X G^,;M^^/3^ 



^/ + Y^(4-^/)^' 



X GIvjM^m/?/^ 



Xl^.xfl^l^f ■ 



487r 

— X Gyj^jM^y^Pf, 

^f X. r^2 A #2 



2 „.2 



47r 



X 



G^jM^m/?/ (7 + Zf) 



where the final state velocity /?/ = ^J\ — Zf and Zf = m'j/M^j^, the color factor Cf 
quarks otherwise 1. The scalar DM and relevant operators av are given by 



aM'ff : 



v/2 ' 



'h^f ■■ 






^ X 2aJ/3j, 
X F^.Ml^Pf 



(2 + ^/) 



^0^^,0/V7V : 



»7r 
An 



SPfl-'t^ 
2 ;\#2 



X F^^fM^^Pf 



(2-^/)' 



(B7) 
(B8) 

(B9) 

(BIO) 
(Bll) 
(B12) 

(B13) 
(B14) 
(B15) 

= 3 for 

(B16) 
(B17) 
(B18) 
(B19) 



Appendix C: A Non-Supersymmetric Scalar iDM Model 

In the non-supersymmetric case, there is a simple way to realize iDM at renormalizable 
level. The dark sector consists of two SM singlets 01^2 with mass hierarchy 0(8 GeV)^ ~ 
m?, <C Tn?j,2- "^^^ relevant scalar potential for iDM generation is quite simple 



-V D KJ0ip + ml\<l)2\') + (r/1010;^^ + V2<l>lS' + h.c) . 



(CI) 



We need to arrange charge assignment properly so that the dark sector conserves the Z2 
symmetry to protect DM stable. Singlet 5* breaks U{l)x at v^ ~ (9(100) GeV as the model 
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given in the text, which also induces a quartic mass term for the heavy state 02 

-V^r,^vl<Pl + c.c., (C2) 

and breaks the dark global ^(1) symmetry (acting on 0j only). It renders small mass 
splitting between the CP-even and CP-odd states of 0i. To extract out the mass splitting 
analytically, we diagonlize the mass matrix for 0i^2 by the unitary matrix t/^, where we 
ignore the U{1) breaking mass term. Then we obtain 

01 -)■ cos 6'i20i + sin ^12025 02 -^ - sin 6'i20i + cos 6'i202, (C3) 

where the fields with prime are in the (approximate) mass eigenstates. And the mixing angle 

^i2=^r/i4-«l (C4) 

is invalid if {rjivlY < mlm^ (assure the positivity) and \riViV2\i'm\ <^ m\. With it, sub- 
stituting Eq. (C3) into Eq. (CI), then it is not difficult to get the mass splitting for the 

01 



6 = 2r/2 sin ^^ x J^ ~ 21^2^1 f ^^1 ^^- (C5) 

Thus, we can obtain 5 ~ 10"^ GeV in many ways, even by setting the m2 at the interesting 
TeV scale. 
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